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ABSTRACT 

Context. The origin of protostellar jets as well as their impact on the regulation of angular momentum and the inner disk physics are 
still crucial open questions in star formation. 

Aims. We aim to test the different proposed ejection processes in T Tauri stars through high-angular resolution observations of 
forbidden-line emission from the inner DG Tauri microjet. 

Methods. We present spectro-imaging observations of the DG Tauri jet obtained with SINFONI/VLT in the lines of [Fe n] /1 1.64 /im, 
1.53/im with 0715 angular resolution and R=3000 spectral resolution. We analyze the morphology and kinematics, derive electronic 
densities and mass-flux rates and discuss the implications for proposed jet launching models. 

Results. (1) We observe an onion-like velocity structure in [Fe n] in the blueshifted jet, similar to that observed in optical lines. 
High- velocity (HV) gas at = -200 km s -1 is collimated inside a half-opening angle of 4° and medium- velocity (MV) gas at ^ -100 
km s~'in a cone with an half-opening angle 14°. (2) Two new axial jet knots are detected in the blue jet, as well as a more distant 
bubble with corresponding counter-bubble. The periodic knot ejection timescale is revised downward to 2.5 yrs. (3) The redshifted jet 
is detected only beyond 077 from the star, yielding revised constraints on the disk surface density. (4) From comparison to [O i] data 
we infer iron depletion of a factor 3 at high velocities and a factor 10 at speeds below -100 km s" 1 . (5) The mass-fluxes in each of 
the medium and high- velocity components of the blueshifted lobe are ^ 1.6 ± 0.8 x 10~ 8 M Q yr 1 , representing 0.02-0.2 of the disk 
accretion rate. 

Conclusions. The medium- velocity conical [Fe n] flow in the DG Tau jet is too fast and too narrow to trace photo-evaporated matter 
from the disk atmosphere. Both its kinematics and collimation cannot be reproduced by the X-wind, nor can the "conical magneto- 
spheric wind". The level of Fe gas phase depletion in the DG Tau medium- velocity component also rules out a stellar wind and a 
cocoon ejected sideways from the high- velocity beam. A quasi-steady centrifugal MHD disk wind ejected over 0.25-1.5 AU and/or 
episodic magnetic tower cavities launched from the disk appear as the most plausible origins for the [Fe n] medium velocity com- 
ponent in the DG Tau jet. The same disk wind model can also account for the properties of the high-velocity [Fe n] flow, although 
alternative origins in magnetospheric and/or stellar winds cannot be excluded for this component. 

Key words. Stars: pre-main sequence - Stars: individual: DG Tau - ISM: jets and outflows - ISM: individual objects: HH158 - 
Techniques: imaging spectroscopy - Techniques: high-angular resolution - 

1. Introduction equipartition magnetic field required to laun ch a powerful M HD 

disk wind could affect planetary migration (Terquem 2003b and 
The physical mechanism by which mass is ejected from accret- wouW speed up ±e accretion flow inside the launch zone low . 

ing young stars and is collimated into jets and the role that eri the disk surface densi b seyeral orderfj of magnitude 

g , these jets play in solving the angular momentum problem of compared to a standard viscous disk model (Combet & Ferreira 

star formation remain fundamental open questions. It is now feppg) 

widely accepted that the action of magnetic fields is r equired 

to exp lain the efficient jet collimation and acceleration dCabritl Strong observational constraints on the jet driving mecha- 
l2007al) . However, several MHD ejection sites may operate in nism in young stars come from studies of extended fordidden 
parallel: the stellar surface, reconnection points in the stellar line emission in "microjets" from young low-mass (< 2M ) 
magnetosphere, the inner disk edge or a broader range of disk T Tauri stars (TTS). Their proximity and lack of obscuring enve- 
radii, and their r elative contribution to the jet remain s subject to lope give access to inner jet regions within 20-200 AU where the 
debate (see e.g. iFerreira et al.ll200"6t lEdwardsl 12009, and refer- flow structure has, hopefully, not yet been significantly perturbed 
ences therein for a recent review). by interaction with the ambient medium. In addition, the opti- 
Distinguishing between these scenarii is crucial not only cally thin forbidden lines provide powerful tracers of jet density 
for understanding the role of jets in star formation, but also and excitation. In particular, a wealth of data on jet widths, kine- 
for planetary formation and disk evolution. For example, the matics, and density and ionization structure has been obtained 
over the last 10 years in the optical range, thanks to the sub- 
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and HST/STIS (see the review bv lRav et al.l2 007. and references 
therein). Confrontin g proposed ejection models with these data, 
iFerreira et al.l d2006l) proposed that centrifugal disk winds ap- 
pear the most promising to reproduce simultaneously all forbid- 
den line jet observations, and that these disk winds could confine 
the hot inner stell ar / magnetqsphe ric winds inferred, e.g., from 
He I line profiles dKwan et al.ll2007l) . 

Support for an extended MHD disk wind is most evident in 
the bright and well-studie d jet HH 158 from the actively accret- 
ing T Tauri star DG Tauri dMundt & Friedll 19831) . An onion-like 
structure is observed within 500 AU of the star in optical for- 
bidde n emission lines, with faster gas nested inside slowe r ma- 
terial (lLavallev-Fouquet et alll2000t iBacciotti et~ai1l2000i) . This 
behavior is naturally expected if the wind is launched from a 
broad range of disk radii. Furthermore, the slowest jet material at 
V r - 30- 60 km s~' exhibits velocity shifts between the two sides 
of the jet (IBacciotti et al.ll2000t ICoffev et al]l2007[) in th e same 
sense as the rotation of the DG Tau disk (Testi et al. 2002). These 
transverse shifts excellently agree with predicted rotation signa- 
tures f or a steady-state MHD disk wind launched out to about 
3 AU dAnderson et al.ll2003t iPesenti et"aTI l2004t IFerreira et all 
2006). A velocity shift in the same sense is seen across faster 
jet material at -200 km s _! suggesting (if caused by rotation) 
ejecti on from smaller disk radii of ^ 0.2-0.5 AU dCoffev et al.l 
120071) . 

However, rotation signatures alone are insufficient to firmly 
establish a disk wind origin for the DG Tau jet, because velocity 
shifts might be contaminated by ot her effects, such as jet preces- 
sion o r non- axisymmetric shocks (Ceraueira et al. 2 006t ISokeil 
120051) . Alternative explanations for the broader, slow component 
in forbidden lines could be, e.g., a cocoon of gas ejected side- 
ways from the central fast jet beam or ambient swept-up gas. To 
further test these alternatives, more information is needed such 
as the full flow morphology in 2D, accurate mass-fluxes, and 
dust content. 

Another open question is whether the fast and slow material 
seen in forbidden lines in the DG Tauri jet share a similar origin 
or if they trace two distinct ejection sites and ejection mecha- 
nisms. A dual or igin in a magnet ospheric wind and a disk wind 
was proposed bv lPvo et alJ ([2003) for the emission in [Fe n]from 
long-slit observations at high spectral resolution (30 km s~') 
based on the double-peaked profile with a high-velocity com- 
ponent at about -200 km s _1 and a "low-velocity" component at 
about -100 km s , of relative strength varying along the jet. In 
contrast, a single disk wind where the ionization fraction varies 
continuously from inner to outer streamlines wo uld predict a 
smoot h single-peaked profile, as shown in Fig. 2 of lPesenti et al.l 
(2004). However, the ionization distribution in the disk wind has 
a crucial effect on the line profile shape, with, e.g., the peak shift- 
ing from high to low velocity d epending on whethe r ionization 
decreases or increases outward (Pesenti et al. 2004j). Therefore, 
double-peaked profiles could in principle be produced in a sin- 
gle disk wind if ionization conditions are highly inhomogenous 
between inner and outer streamlines, e. g. because of localized 
shocks. To better constrain the respecti ve origins o f the t wo ve- 
locity components noted in [Fe n] by iPvo et al] d2003l) . addi- 
tional constraints are again needed. 

In this paper, we present spectro-imaging of the DG Tau jet 
in the H-band that provide the first images of this jet in [Fe n] at 
— 0'.' 15 resolution over a full field of view of + l'/5 from the star. 
This enables us to clarify the morphology and openi ng angles 
of the two [Fe n] velocity components identified by IPvo et alJ 
(2003) in long-slit spectroscopy. We also identify new knots 
and bubble features in both the jet and counterjet, and derive 
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Fig. 1. Illustration of the continuum substraction methods, close 
to the star (bottom row) and far from the star (top row), for the 
[Fe n] 1 .64/im (left panels) and the H2 2. 12yum line (right panels). 
The raw spectra are plotted in thick lines and the fitted continua 
in thin lines. See text for details. The analysis and discussion of 
H2 2. 12yum lines will be presented in a separate paper. 

the electronic density distribution along the jet, using the [Fe 11] 
1.53/^m/1.64//m ratio. The observations and data reduction are 
presented in Section [2] and our results in terms of morphol- 
ogy, kinematics, collimation, red/blue asymmetry, disk occulta- 
tion, and jet densities are presented in Section [3] In Section |4] 
we use these results to revise the timescale for knot formation 
in the DG Tau jet, measu re iron depletion in the two velocity 
components identified by IPvo et al.l d2003) and estimate their 
mass-flux and ejection to accretion ratio. These results are com- 
pared to proposed ejectio n processes in T T auri stars, in par- 
ticular the X-wind model dShang et alJI 1 998b and the extended 
MHD dis k wind model invo ked to fit the rotation signatures in 
DG Tau (Pese ntiet al.1 12004). Our conclusions are summarized 
in Section[5] 

2. Observations and data reduction 

Spectro-imaging observations of the DG Tau microjet in the 
H- and K-bands were conducted on 2005 October 15 at the 
Very Large Telescope (VLT) with the integral field spectro- 
grap h SINFONI combined with an adaptive o ptics (AO) mod- 
ule (lEisenhauer et al.ll2003UBonnet et al.ll2004l) . We used a spa- 
tial sampling of C/05 xO'.'l, with the smallest spaxel dimen- 
sion aligned at PA = 315°, i.e. in the direction transverse to 
the jet axis, denoted x-axis in the following. This provided a 
total field of view of about 3" x 3". The spectral resolution 
was AA/A -3000 in the H-band and AA/A -4000 in the K-band, 
with a spectral sampling of 35 km s _I /pixel and 34 km s^ 1 /pixel 
in H and K respectively. In the H-band, individual short expo- 
sures of 15sec were coadded to build a final datacube with to- 
tal integration time of 2x270s on the object and 270s on the 
sky. Thanks to AO correction, the effective spatial resolution 
achieved is ~0.15"FWHM in the PSF core (estimated from the 
reconstructed continuum image). 

The data reduction process is common to the two spectral 
bands. Data reduction steps were carried out with the SINFONI 
pipeline. Datacubes were corrected for bad pixels, dark, Hatfield, 
geometric distortions on the detector, and sky background. The 
combined datacube was resampled to a 50 mas x 50 mas square 
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spatial grid. A raw wavelength calibration based on daytime ar- 
clamp exposures is performed by the pipeline. However, the ac- 
curacy of this calibration (on the order of 50km s -1 ) is insuf- 
ficient for our kinematical jet studies. For the [Fe n] 1.64yum 
line, we improved this wavelength calibration locally using the 
bright OH sky emission line at /ll.6692/im to obtain a 2D wave- 
length solution over our whole field of view. We estimate a final 
calibration error in velocity of 10km s _1 rms in this line. The 
wavelength scale was converted into the stellar rest frame veloc- 
ity scale. Unfortunately, for the [Fe n]/ll.53;um line, nearby OH 
lines have a signal-to-noise ratio that is too low to improve wave- 
length calibration. Therefore, only velocity-integrated 1.53/im 
line fluxes are used in the analysis. The data were flux-calibrated 
on the A2V standard star HIP 26686 in the H-band. We esti- 
mate an error in the flux calibration of 5%. Correction for tel- 
luric absorption lines was performed together with our contin- 
uum subtraction routine (see below). To correct for differential 
atmospheric diffraction and ensure proper registering, the final 
datacube for each line of interest is spatially recentered on the 
centroid position of the local continuum. 

The bright stellar spectrum is strong over most of the field 
of view. In particular, a strong photospheric line is located un- 
der the [Fe n] 1 .64/im transition. To better retrieve the intrinsic 
[Fe n] jet emission profile, especially close to the central source, 
we performed a specific continuum subtraction, illustrated in the 
bottom row of Figure [1] We selected a high signal-to-noise refer- 
ence photospheric spectrum free of jet emission in the SINFONI 
field of view (at Ax=0.054", Ay=-0.003"). This reference spec- 
trum was then scaled down to fit the local continuum level and 
was subtracted out from the current spectrum. The high qual- 
ity of this continuum-fitting procedure can be seen in the bottom 
row of FigureQ] It allows us to retrieve the extended jet emission 
features down to 072 from the star. Furthermore, the residual 
emission line profile is not only free of photospheric absorption 
lines but also of telluric features against the continuum, because 
the depth of telluric absorption scales in proportion to the local 
flux leveQ Farther from the star, where the line-to-continuum 
ratio becomes high, this procedure only provides a partial cor- 
rection for telluric absorption. However, the telluric absorption 
depth is less than 7% locally around the [Fe n] lines of inter- 
est, therefore the error in the emission-line flux is comparable 
to the flux-calibration uncertainties. At distant positions, where 
the signal to noise in the stellar continuum is too low for proper 
continuum scaling, a simple linear baseline is fitted to the local 
continuum around the line and subtracted (top row in Figure [TJ. 
At each position, we also estimate a spectral noise from the stan- 
dard deviation of the residual (continuum-subtracted) spectrum 
on either side of the line. The same continuum subtraction pro- 
cedures were applied to the H2 2.12/im line in the K-band. The 
results in H2 will be presented in a separate paper. 

3. Results 

The present H-band SINFONI data provide the first observations 
on sub-arcsecond scales of the 2D spatio-kinematic structure of 
[Fe 11] emission in the DG Tau jet. Below, we present in turn 
several new results obtained from our data: 

- the [Fe n] jet morphology and collimation in various velocity 
intervals 

1 Our photospheric subtraction procedure would also subtract any 
unresolved [Fe 11] feature centered on the star. However, no such com- 
pact co mponent appears to be present in previous high-resolution [Fe 11] 
spectra dPvo et alj2003h . 



- the kinematic structure and variability of [Fe 11] emission, 
compared to earlier observations in [Fe n] and optical lines 

- the red/blue asymmetry and constraints on the occulting disk 

- the electronic density distribution derived from the [Fe 11] 
1 .53/im/l .64yum ratio. 

3.1. 2D morphology in [Fe 11/ 

We show in the top row of Fig.[2]the raw continuum-subtracted 
[Fe 11] channel maps reconstructed from our datacube. The emis- 
sion in each jet lobe was separated in two velocity ranges, cho- 
sen to correspond best to the two kinematical comp onents at - 
200 km s" 1 and -100 km s _1 spectrally identified by iPvo et al.l 
(2003) in their 2001 [Fe 11] long-slit spectrum along the jet 
axis Q. The velocity intervals are the following: high-velocity 
blueshifted (HVB) = [-300,-160] km s -1 , medium-velocity 
blueshifted (MVB) = [-160,-50] km s _1 , medium-velocity red- 
shifted (MVR) = [-50,120] km s _1 , high-velocity redshifted 
(HVR) = [120,260] km s _1 . Because of our moderate velocity 
resolution and the lower velocities in the redshifted jet, we ex- 
tended the MVR channel map to -50 km s _1 . The HV channel 
map shows the sum of the HVB and HVR intervals. The MV 
channel map shows the sum of MVB and MVR. The blueshifted 
lobe is at Ay < in Fig. [2] 

The bottom row of Fig. [2] shows the channel maps after 
deconvolution by a continuum image reconstructed from the 
same datacube, which provides a good estimate of the point- 
spread function. We used the LUCY algorithm implemented in 
the IRAF/STSDAS V3.8 with 40 iterations. We checked that 
the derived jet morphology (in particular the jet widths) did not 
change significantly with further iterations. The right-hand panel 
of Fig. |2]presents a color superposition of the deconvolved HV 
and MV channels to facilitate comparison between them and to 
highlight faint emission features. 

The deconvolved channel maps reveal several important new 
features: 

(1) the blueshifted [Fe 11] jet shows a clear increase in colli- 
mation at higher flow velocities, with the HVB emission nested 
inside the MVB emission. This on ion-like velocity structure was 
previously noted in optical lines (lLavallev-Fouquet et al.l l2000: 
iBacciotti et ali 2000), but the unique combination of high resolu- 
tion and wide transverse coverage provided by SINFONI shows 
it more clearly. In particular, it confirms the striking conical 
shape of the MVB emission, which was partly truncated in the 
narro wer HST/STIS chan nel maps covering +072 from the jet 
axis (IBacciotti et al.ll2000l) . 

(2) The blueshifted jet contains several knots and bubble-like 
structures: Two knots are present, one at -0737 +0.03" and one 
at -1'.'2 ±0.05". The latter is followed by a faint limb-brightened 
bubble, partly truncated by our field of view, which reaches a 
maximum width 075 around d -174 (see right panel in 
Fig. [2). 

(3) The counterjet emission is dominated by a striking limb- 
brightened circular bubble centered around +1". The decon- 
volved MVR emission is slightly broader and more pointed away 
from the star than the HVR. No counterjet is detected within 077 
from the star. This effect is probably c aused by disk occultation 
(lLavallev et al.|[l997t IPvo et al] 12003) and provides interesting 



2 We will use here the notation MV ( medium- ve l ocity ) rather than 
LV for the component at =s -100km s -1 in lPvo et all (2003) to be more 
co nsistent with the vel ocity intervals defined in the HST channel maps 
of IBacciotti et al. (2000); we restrict the use of LV to the component at 
velocities lower than -50 km s _I , seen in optical lines. 
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Fig. 2. Left : Continuum-subtracted [Fe n] 1.64/im channel maps of the DG Tau jet before (top panels) and after (bottom pan- 
els) deconvolution by the continuum image. The HV channel maps are integrated over the velocity intervals [-300,-160] km s _1 
(HVB) and [+120,+260] km s 1 (HVR); the MV channel map over [-160,-50] km s" 1 (MVB) and f-50,+120] km s 1 (MVR). The 
blueshifted lobe is at Ay < 0. The pixel scale is 0705. A cross denotes the centroid of the continuum emission (shown as an insert 
in the top left panel). Contours start at 3 cr and increase by factors of V2 (lowest contours are 2.0 and 1.8 xl0~ 18 W irT 2 arcsec" 2 
in the raw HV and MV maps respectively). The solid red lines in the deconvolved maps show the opening angle measured in each 
component. Right : Superposition of the deconvolved HV (background color image) and MV (yellow contours) channel maps. The 
stellar continuum is shown in red contours. The white cross locates the continuum centroid. 



constraints on the disk surface density structure in DG Tau (see 
SectionEy}. 

3.2. Collimation, opening angle, and upper limits on launch 
radii 

We plot in the left panel of Fig. [3] the jet FWHM derived from 
Gaussian fitting of the transverse intensity profiles in the decon- 
volved channel maps of Fig. [2] The broader width of the MV 
channels compared to the HV is clearly visible at all distances. 
In addition, a widening of the blue jet is seen beyond 1", corre- 
sponding to the faint bubble feature at the end of this lobe. The 
width there is comparable to the bubble feature in the redshifted 
lobe. 

To estimate intrinsic jet widths and opening angles, we fo- 
cused on the inner 078 (110 AU) section of the blueshifted 
lobe, which is least affected by bubble structures. We corrected 
the measured FWHM for the effective PSF width (FWHMo) 
in the deconvolved channel maps with the formula (2Rj et ) 2 = 



FWHM 2 - FWHM 2 ,. The effective angular resolution reached 
in the deconvolved maps cannot be measured precisely be- 
cause no unresolved point source is present in our (continuum- 
subtracted) line maps. Therefore we considered a range of pos- 
sible FWHM between 0705 (pixel scale) and 0708 (smallest jet 
width measured). We plot in the right panel of Fig. [3] the PSF- 
corrected jet radii (in AU) in the HVB and MVB channels as a 
function of distance as well as linear fits to these data points. 

We find that the medium-velocity blueshifted (MVB) emis- 
sion has a conical shape with a roughly constant half opening an- 
gle of 14°. The high-velocity blueshifted emission (HVB) keeps 
a narrow radius (5-10 AU) with a half opening angle of 4° until 
0.8" from the star. The derived opening angles are upper limits 
because we did not correct for inclination effects. If these cones 
traced actual streamlines, an extrapolation back to the origin 
would suggest that the HVB component in DG Tau originates 
in a disk radius < 4 AU and the MVB in a disk radius < 1 AU. 
This is compatible with the maximum launch radii of 0.5 AU 
and 2-3 AU derived from tentative jet rotation signatures at pro- 



V. Agra-Amboage et al.: Sub-arcsecond [Fe n] spectro-imaging of the DG Tau jet 



5 



> MVR ,-'\ 




□ 



DG Tau [Fe II] HVB 
DG Tau [Fe II] MVB 

HN Tau, UZ Tau E - 
Hartigan et al (2004) 

RW Aur - [0 |],[S 111 HVR 
Woitas et al (2003) 

DG Tau - [0 I] All Velociti 
Dougados et al. (2000) 




Projected Distance to the star (in arcsec) 



Projected Distance to the star (in AU) 



Fig. 3. Left: Jet FWHM derived from Gaussian fits to deconvolved [Fe n] channel maps as a function of projected angular distance 
from the star. Right: PSF corrected (i.e. intrinsic) jet radii as a function of linear distance from the star for the inner 078 of the HVB 
(solid curves) and MVB (dashed curves). The fitted FWHM across the jet were corrected for the angular resolution of 0705-0708 in 
the deconvolved maps by subtracting this PSF in quadrature. Gray/green solid lines show fits by a flow of constant opening angle, 
suggesting a launch radius smaller than 4 AU (HVB) and 1 AU ( MVB). For comparison , open symbols plot j et radii measured i n 
[O i] images integrated over all velocities in DG Tau (open squar es Dougado s et al.ll2000h . RW Aur (open stars IWoitas et al.l l2002). 
and HN and UZ Tau E (open circles and triangles Hartigan et al. 20 04al) . 



jected flow speeds of -200 km s 1 and — -50 km s 1 respec - 
tively (ICoffev et al.l l2007h iBacciottil l2002t iPesentiet all 120041) . 
However, it is likely that emission maps in constant velocity bins 
do not trace actual streamlines especially at low flow velocities. 
Indeed, outer streamlines are expected to reach their terminal 
velocities on spatial scales ^ 50-100 AU. As a consequence, jet 
emission radii measured on the MVB channel maps likely un- 
derestimate true streamline radii close to the source. 

We also compare in Fig. [3] our derived intrinsic jet radii 
in [Fe n] with previous measurements from velocity-integrated 
[O i] images at similar angular resolution in the RW Aur jet 
(average of both lobes) (IWoitas et alj|2002|) and in th e HN Tau 
and UZ Tau E blueshifted iets (iHartigan et al.ll2004ah . The RW 
Aur microje t, where \Q il is dominated by a single high- velocity 
component ( Woita s et al.l 12002b . shows an opening angle and 
width strikingly similar to the HVB interval in DG Tau. The 
HN Tau and UZ Tau E [ O il jets, dominated by a medium 
to low-velocity component (Hartigan et al. 1995), have a much 
wider opening angle similar to the MVB. Therefore, the increase 
in collimation at higher velocity seen in DG Tau could be a 
widespread generic property in T Tauri jets. Finally, we note that 
the jet widths from a ve locity-integrated \Q il image of the DG 
Tau jet taken in 1997 dDougados et all 120001) fall close to the 
MVB component out to 50 AU, but closer to the HVB beyond 
50 AU. Therefore, jet opening angles measured from velocity- 
integrated images may be misleading when different velocity 
components dominate at different positions along the jet. 

3.3. 2D kinematics of the blueshifted jet 

To attain an overall view of the 2D jet kinematics, we performed 
at each spatial position a one-component Gaussian fit to the 
line profile. Although two kinematical components are clearly 
present in the jet (as shown by the distinct morphologies of chan- 
nel maps in the HV and MV velocity ranges), a two-component 
Gaussian fit did not give coherent results over the field of view 
because of the moderate spectral resolution of our observations. 
We show in Fig. |4] the derived 2D map of the fitted centroid 



veloc ities. From the formula outlined in lAgra- Amboage et al.l 
(2009) and the observed signal-to-noise ratios, we estimate an 
uncertainty in centroid velocities in the range 14-21 km s _1 , 

The centroid map shows a clear transverse velocity gradi- 
ent in both lobes, with the highest flow velocities concentrated 
toward the axis. The gradient is especially apparent in the blue 
lobe at distances greater than 0.4" from the star. This transverse 
gradient is consistent with the higher degree of collimation in the 
HV channel maps and the wider, limb-brightened morphology 
in the MV channel maps. The centroid values toward the flow 
edges indicates a mean absolute velocity of < 50 km s _1 for the 
MVR component, and < 100 - 150 km s _1 for the MVB com- 
ponent toward the blue l obe, similar to th e low velocity spectral 
component identified by IPvo et al.l (120031) in their 2001 long-slit 
[Fe n] spectrum, of higher spectral resolution than our data. 

Variations of centroid velocities along the jet axis are also 
seen. They are illustrated in Fig. [4] right panel, where we plot a 
position-velocity diagram along the jet, reconstructed by averag- 
ing our spectra over a 1" wide pseudo-slit. Transverse-average 
centroid velocities decrease from -192 km s~'at -0.4" (posi- 
tion of the first knot) down to -150km s~'at -0.8" and increase 
again to reach -210 km s _1 at the outer knot at -1.2". We note 
a striking similarity with the amplitude and period of the veloc- 
ity var iations in the high-velocity range reported by IPvo et all 
(120031) along the jet, despite the time difference of four years 
between our observations. When corrected for a jet inclination 
of 42°, the variation amplitude AV ra d ^ 40-60 km s sug- 
gests internal shock speeds v s ^ 60 - 80km s _1 , similar to the 
v s ^ 50 - 70km s _1 inferred in 1998 from optical fine ratios 
in the velocity range [-250,-100] km s _1 (Lavallev-Fouauetet al. 
2000). In both lobes, the maximum radial velocity is found to- 
ward the bubble feature near the edge of our field of view. 

One notable d i ffere nce between our [Fe n] PV diagram and 
that of IPvo et alJ ([2003) is that the line flux is dominated by 
the HVB all the way down to -072 from the star in our obser- 
vations (Fig. |4] right panel), wherea s it was domina ted by the 
MVB within 075 of the star in 2001 (IPvo et al]l2003l This dif- 
ference could be caused by the bright high-velocity knot present 
at -074 in our SINFONI map, which had not yet — or barely — 
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Fig. 4. Left: 2D maps of the [Fe n] 1.64/mi velocity centroid deduced from one component Gaussian fitting of the line profile. 
Redshifted velocities are shown on a negative scale to better highlight the velocity asymmetry between them. Spectra with SNR 
lower than 5 were masked. The cross marks the location of the continuum. Contours show the sum of the deconvolved HV and MV 
channel maps. Right: PV diagram of [Fe n] along the jet, averaged over ±075 across the jet, with velocity centroids indicated as 
crosses. Horizontal dash-dotted lines indicate the region where residual noise after continuum subtraction is still large. Contours 
start at 2.43xl0~ 15 W rrT 2 /irrr 1 arcsec~ 2 and increase by factors of y2. The 3cr contour is shown as a dotted line. 



emerged in 2001 October (see Section 147Tb . Another difference 
is the nondetection of th e MVB component beyond 075 from the 
star bv lPvo et al.l J2003). which may be caused by the narrow slit 
(073 width) used in their observations. As shown in Figure [3] 
the measured MVB FHWM exceeds 073 beyond 076, therefore 
a significant fraction of the MVB flux has likely been missed at 
these distances in the lPvo et alj d2003l) long-slit observations. 

Finally, we note that neither our PV diagram nor that of Pyo 
et al. from 2001 show [Fe n] emission at velocities faster than - 
300 km s~'or lower than -60 km s , corresponding respectively 
to the "very high velocity" (VHV) and "low velocity" (LV) com- 
ponents dete cted in the DG Tau blueshifted jet in optical lines 
respectively (|Lavallev et al1ll997tlLavallev-F ouquet et al. 2000; 
iBacciotti et al.ll2000l) . HST maps obtained in 1999 show that the 
VHV optical gas was tracing a bri ght knot followed by an elon- 
gated bubble (Bac ciotti et al.ll2000l) . It was therefore produced 
by an intrinsically transient event, and it is conceivable that new 
knots probed by the 2001-2005 [Fe n] observations may not 
have reached these high velocities. 

On the other hand, the halo of broad LV emission seen in 
[O i] and [S n] surrounding the blueshifted DG Ta u jet is present 
in all optical observations since at least 1993 dLavallev et al.l 
119971: iLavallev-Fouauet et al l 12000] IBacciotti et al.l l2000h and 
most recently in 2003 dCoffev et al.ll2007l) . Therefore it appears 
to be an intrinsically stable feature, and its absence in [Fe 11] is 
significant. To illustrate this, we compare in Figure |5]the trans- 
verse PV diagrams of [S 11] 6731 A and [O 1] 6300A obtained 
on 2003 December 1 with HST/STIS at -073 from the source 
dCoffev et al.l l2007) with a transverse PV diagram in [Fe 11] cov- 
ering the same area, reconstructed from our SINFONI datacube. 
In optical data, a spatially broad LV component is clearly visible 
of the centroid decreasing from -6 km s~ l on-axis t o -30,-40 
km s" 1 at ±072 from the jet axis dCoffev et al.ll2007h . Similar 
low speeds were observed at least out to 076 from the star (see 
parallel slits SI and S7 in Fig. 3 of BacciottH l2002l) . The weak- 
ness/absence of [Fe 11] emission from this "off-axis" LV com- 



ponent is striking. In Section [4] we interpret this difference as 
substantial iron depletion onto dust grains in the low-velocity 
sheath surrounding the jet. 



3.4. The counterflow: velocity asymmetry and constraint on 
disk structure 

The counterjet emission in our data is dominated by a striking 
bubble-like structure of similar width and brightness as the bub- 
ble partially imaged in the blueshifted jet (see Figs. 1213b . On the 
other hand, the velocity centroid map presented in Fig.|4]reveals 
a clear velocity asymmetry between the red and blue bubbles. 
The maximum measured radial velocity centroid toward the red 
bubble is ^ +160 km s _1 , ie. only 70% of the peak velocity ^ 
-230 km s _1 in the blue bubble. This velocity asymmetry is con- 
sistent with the 70% smaller distance from the star to the center 
of the redshifted bubble, compared to the (approximate) center 
of its blueshifted counterpart (^ +1" versus -174). This agree- 
ment suggests that both bubbles were generated in a simultane- 
ous ejection event. 

The red/blue asymmetry in velocity and position needs 
to be taken into account when quantifying the dust extinc- 
tion produced by the DG Tau disk toward the counterjet. In 
October 2001. iPvo et all d2003l) detected only very wealfl red- 
shifted [Fe 11] emission around 079, 10 times weaker than their 
blueshifted jet at the same distance from the star. They inferred 
an extinction of 2.5 mag at H at this distance, assuming in- 
trinsically equal brightnesses. If one corrects for the differen- 
tial speeds and compares instead with their blueshifted jet inten- 
sity at 079/0.7 = 173, the relative intensities become comparable, 
suggesting negligible extinction by the disk at 1". This result is 



3 the redshifted bubble was probably at that time entirely hidden be- 
hind the high-extinction region close to the star. Assuming a proper 
motion of 70% that of the blueshifted knots (^ 073/yr, see Section RTH 
the tip of the redshifted bubble projected at 1'.'3 arcsec from the star in 
October 2005 would have been at 0.'46 in October 2001. 
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Fig. 5. Comparison of transverse PV diagrams of 
[Sn]673lA and[Oi]6300A obtained in 20 03 with HST/STIS 
atd=-073 from the source in a 071 wide slit dCoffev et al.ll2007L 
2008) with a transverse PV diagram in [Fe n] at the exact same 
position, reconstructed from our SINFONI datacube obtained 
in Oct. 2005. Note the striking lack of [Fe n] emission in the 
spatially wide low-velocity [O i] component at V=-30,-40 
km s . 



much more consistent with the similar brightness of the red and 
blue bubbles in our SINFONI data. 

The deconvolved [Fe n] channel maps suggest that the cen- 
tral obscuration of the redshifted jet from the disk extends to 
projected angular distances of d ^ 0.7", i.e. 140 AU in depro- 
jected disk radius (assuming a distance of 140 pc and a disc in- 
clination of 42°). This requirement appears to set constraints on 
models of t he dis c density distribution in DG Tau derived by 
Isella et al.l d2010t) from sub-arcsecond millimeter continuum in- 
terferometric maps. In particular, it seems inconsistent with a 
single power-law surface density distribution, which requires a 
small outer disk r adius of 70-80 AU to fit DG Tau observations 
(llsella et al.ll2010h . Such a small disk radius would produce no 
extinction beyond 074 from the star, unlike SINFONI observa- 
tions. An alternative similarity solution with an exponential tail 
fitted to the mm map predicts a (perpendicular) surf ace density 
of 2=0.25-1.5 g cirT 2 at a disk radius of r=100 AU (llsella et all 
l20loh . leading to Ah = 30-5 mag at d=0.5" assuming interstel- 
lar dust extinction, with X and therefore Ah dropping by a factor 
s 10 at d=0.7" (r=140 AU). The similarity solution for the sur- 
face density distribution therefore appears more consistent with 
the central obscuration observed in the SINFONI maps. 

3.5. Electron density in the [Fen] emitting region 

The ratio R of the [Fe n] 1.53yum to 1.64/im lines is a good 
diagnostic of the electronic densit y n e in the range 10 3 - 
10 5 cm -3 dPradhan & Zhang! Il993l) . The temperature depen- 
dence of R is weak in the range T e =s 8000 - 20,000 K encoun- 
tered in the forbidden-line emission regions of stellar jets (e.g. 
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Fig. 6. Plot of n e versus distance derived from the ratio of the 
[Fe n] 1.53/mi to 1.64/mi lines (averaged over +075 from the 
jet axis) for the blueshifted (blue points, top curve) and red- 
shifted (red points, bottom curve) lobe. Also shown are previous 
determinations from the [S n] d oublet ratio over similar veloc- 
ity ranges: IVB and LVB from lLavalle v-Fo uq uet et al. (200< 

black circle, purple square), MVB interval of iBacciotti et al 

2000) (green triangles). 



Bacciotti & Eisloffel 119991) . We did not detect any change in this 
ratio across the line profile nor across the jet width at our reso- 
lution. Therefore we computed at each spatial distance along the 
jet a mean [Fe n] line ratio integrated over all velocities and over 
±075 from the jet axis. The flux for each line was obtained by 
fitting it to a Gaussian and then integrating over all velocities. 
Errors in the Gaussian parameters are derived taking into ac- 
count observational errors on the spectral profile and propagated 
to obtain the final errors on the derived ratio. We then derived the 
electronic density from the observed line ratio using the theoret- 
ical relation between R and n e at a temper ature T e ^ 10 , 000 K , 
calculated from the 16-level Fe + model of iPesenti et al . (2003). 
For ratios R between 1.0 and 3.9, this relation can be fitted ana- 
lytically within an accuracy of 20% on « e by 



n e (R) = 1 .2 X 10 4 X (R - 0.035) / (0.395 - R) cm" 



(1) 



The n e values as a function of distance along the jet are pre- 
sented in Fig. [6] For comparison, we also plot in Fig. [6] earlier 
« e determinations obtained in January 1998 and 1999 from the 
optical [S n] 67 16,673 1AA doublet ratio over a similar velocity 
range: namely, the IVB [-250,-100] km s"' and L VB [-100,+ 10] 
km s" 1 intervals o flLavallev-Fouauet etall (l2000t) and the MVB 
[-195,-72] km s" 1 interval of lBacciotti et al 1(12000) . The HVB [- 
319,-195]km s _1 interval of Bacciotti et al. had a saturated ratio 
indicating electronic densities » 10 4 cirT 3 out to 174. 

This comparison shows that the overall n e distribution along 
the blueshifted DG Tau jet has remained remarkably similar in 
the velocity range [-300,+10] km s _1 between 1998 and 2005, 
with a general increase of only a factor 2 over this period. The 
actual change in total flow density, n#, would be even smaller 
than a factor 2 if the ionization fraction had increased over this 
time period or if the [S n] 671 6,673 1AA doublet were biased 
toward less dense jet regions than [Fe n], as found bv Nisin i et al .1 
(2005). We will take advantage of this small variability level to 
estimate iron gas-phase depletion in Section l4~2l 
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the D G Tau blueshifted jet updated from Fig. 5 of 
( 2003). Solid lines show the proper motions derived by 
d2003l) for knots Al and A3. For knots with only one epoch, 
dashed lines show the typical proper motion range of 0/27 yr~'- 
0/34 yr . The two new knots observed with SINFONI and 
knot A2 observed with HST indicate that a new knot typically 
emerges every 2.5 years since 19 91. References for knot posi- 
tions may be found in Table 1 of Pvo et aTj (120031) . 



4. Discussion 

Below, we analyze our SINFONI results and discuss the result- 
ing implications concerning 

- the timescale for knot and bubble formation 

- the gas-phase iron depletion in the blueshifted jet 

- the jet mass-flux and ejection to accretion ratio 

- the proposed ejection processes and launch radii in DG Tau. 

4.1. Period of knot generation in the DG Tau jet 



IPvo et all d2003l) investigated the knot proper motions and ejec- 
tion timescales in the DG Tau blueshifted lobe from 1991 to 
2001, using a compilation of data from the literature. They con- 
cluded that a new knot emerges every five years. Here we re- 
examine this matter and show that the SINFONI data indicate 
half that timescale. 

The HVB emission map reveals two well defined [Fe n] 
emission knots at distances from the star of -0737 + 0.03"and 
— 1'.'2 ± 0705. In Figure [7] we add these two HVB knots onto 
the distance-date diagram of IPvo e t al. (2003) (th eir Fig. 5). We 
also plot the proper motions determined by Pvo et al.l d2003l) for 
knot Al (07275 yr -1 from four epochs) and for their HV peak 
(0734 yr~' from two epochs). We denote the lat ter A3, extend- 
ing th e knot-naming convention introduced by Bac ciotti et al.l 
d2000h . 

Evidently, the two high-velocity knots identified in our HVB 
map correspond to new ejection events: the HV knot A3 should 
lie at =s -2" at the time of our observations, i.e. well outside of 
the SINFONI field of view. If we assume for these new knots 
A4 and A5, and for the previous knot A2, proper motions in the 
same range 0727-0734 yr~' as found for Al and A3, we derive 
intervals between knots of about 2.5yrs between 1995 and 2005 
(typical spacing between knots of 0.75 " in the blue lobe). Our 
observations therefore suggest the presence of a knot variabilit y 
timescale a factor two shorter than derived bv lPvo et al.l (|2003). 



We note that the bubble structures observed in our SINFONI 
images appear to have formed in the wake of the same ejection 
event as HV knot A3: if we adopt a proper motion in the red 
lobe of 70% of the mean proper motion of 073 yr~' in the blue 
lobe, as suggested by the red/blue asymmetry (see Section l3~4"l >. 
the position of the tip of the red bubble at +173 from the star in 
October 2005 yields an ejection date around August 1999, the 
same as for knot A3. The formation of a bubble structure in the 
wake of a knot has been previously s een in the case of knot Al 
in HST images dBacciotti et al.ll2000h . It may thus be a common 
feature in DG Tau. 

We did not include the MV knot identified by IPvo et all 
(2003) in our timescale analysis because the estimated proper 
motion of 0728 yr -1 yields a ratio of proper motion to radial 
velocity 60% higher compared to the HV knot. This is not con- 
sistent with ballistic motions of a gas parcel along the jet axis, 
which would predict a constant ratio of proper motion to radial 
velocity (equal to tan(/), with i the jet inclination to the line of 
sight). Therefore the MV knot of Pyo et al. does not appear to 
trace a separate ejection event, but rather a brightness enhance- 
ment following the HV knot A3. 



4.2. Iron depletion in the blueshifted jet 

The dust content of stellar jets is a key parameter in determin- 
ing whether the outflowing material originates inside or outside 
the dust sublimation radius. The gas phase depletion of refrac- 
tory species such as Fe or Ca, measured with respect to the so- 
lar abundance, offers a direct indication of t he presence of dust 
in the jet. In previou s studi es of HH jets bv INisini et al.l d2002i 
l2005l) : |PesentietalJ d2004l) : lPodio et al.l d2006ll2009l) . relatively 
high Fe and Ca gas-phase abundances of 20%-80% the solar 
value were inferred, but most of these observations probe rela- 
tively large distances from the driving sources > 10", where the 
original dust content of the flow may have been m odified, eg. 
throug h du st destructio n in sh ocks. Exceptions are INisini et al.l 
(2005j) and iPodio etail d2009l) who reported Fe and Ca deple- 
tions at 1-3" of the source in the HH 1 and HH 24 jets, cor- 
responding to 450-1300 AU at the distance of Orion. Our sub- 
arcsecond resolution SINFONI observations of the DG Tau jet in 
[Fe n] provide the first opportunity to constrain the dust content 
of the blueshifted jet down to only 50 AU (073) from its driv- 
ing source, i.e. much closer than achieved in previous depletion 
studies of HH jets. 

Iron depletion is estimated by comparing [Fe n] fluxes with 
lines from spe cies that show little or no depletion onto grains , 
such as s ulfur dBeck-Winchatz et al. l l 19961: iPesenti et aljJ2003l). 
oxyge n dMouri & Taniguchil 120001). hydrogen dNisini etldl 
2002), or phosphorus dPodio et al.l2006l) . Because our SINFONI 
data include no lines from the above species, we use here pre- 
vious flux-calibrated [O i]/16300A observations of the inner re- 
gio ns of the DG Tauri blueshifted jet obtained in January 1998 
by lLavallev-Fouquet et al.l d2000h with OASIS/CFHT with an 
angular resolution of 074 and a velocity resolution of 90 km s . 



We show in Fig. |S] the [Fe n] 1.64//m/[0 i] 6300A line ra- 
tios derived from integration of our data over the same veloc- 
ity intervals IVB[-250,-100]km s~'and LVB[-100,+10]km s _1 
and the same distance ran ges along the jet as used by 
lLavallev-Fouauetetal.1 d2000h : d// = -07225, -0/725, -17225, 
with pixel size of 0/5 along the jet. These ratios are plotted 
against the 1 .53/mi/l .64/im ratio, averaged over the same veloc- 
ity and distance intervals. 
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Fig. 8. Line ratio diagnostic diagram for Fe gas phase depletion: 
[Fe n] 1.64/im/[0 i] 6300A versus [Fe n] 1.53/1.64 /an. [Fe n] 
line fluxes are integrated o ver the same velocity and spati al in- 
tervals as [O i] fluxes from lLavallev-Fouquet et al.l d2000). The 
three data points correspond to distances along the DG Tauri 
blueshifted jet of 0725, 0775 and 1725 (increasing with decreas- 
ing [Fe n] 1.53/1.64 fim). Full curves show predictions for vol- 
ume emissi on at 8000 K (top pan el) and for J-shocks (bottom 
panel) from Hartiga n et al.l (|2004bl) assuming an [Fe/O] solar gas 
phase abundance of 0.062 (Asplu nd et al.l|20 05). Dashed curves 
show [Fe/O] abundances required to reproduce observed line ra- 
tios. They indicate Fe gas phase depletion by a factor 3-4 for the 
IVBt^SO^lOOJkms- 1 and 10-12 for the LVB[-100,+10]km s" 1 
intervals. 



In the bottom and top panels of Fig. [8] we compare there 
with predictions for a solar abundance ratio Fe/O = 0.062 
(Aspl und et al.ll2005h in two excitation cases: (i) a stratified pla- 
nar shock wave (denoted "shock method" hereafter), and (ii) an 
emitting volume of homogeneous density and temperature (de- 
noted "volume method" hereafter). 

For at omic shock w a ves, w e considered the models pre- 
sented in lHartigan et alJ (12004b). covering preshock densities 
«h = 10 3 - 10 5 cm 3 , shock speeds of 30-50 km s _1 , a preshock 
ionization fraction of l%-50%, and a preshock magnetic field of 
30-200^/G. We find that the observed [Fe n]/[0 i] line ratios in 
the DG Tau jet are systematically lower than the predictions for 
a solar Fe/O ratio, by a factor 3 in the IVB interval and by a fac- 
tor 10 in the LVB interval, independent of the distance from the 
source. Note that the predicted line ratio curve assumes that the 
mass-flux entering the shock was the same in 1998 ([O i] data) 



and 2005 ([Fe n] data), i.e. that the twice higher electronic den- 
sity in [Fe n] in 2005 compared to [S n] in 1998 (Section [331) 
is mostly caused by ionization gradients within the shock cool- 
ing zone. This gives a conservative lower limit to iron depletion. 
If, instead, the increase in n e between the [O i] and [Fe n] ob- 
servations were caused by an increase in pre-shock density, one 
would predict a higher [Fe n]/[0 i] ratio and infer a larger iron 
depletion in the jet. 

For the uniform volume method, assumptions need to be 
made on the temperature and ionization state of the gas: we as- 
sume a temperature of 8000 K typical of stellar jets, that iron is 
all in the form of Fe + and oxygen is mostly neutral (as expected 
at such temperatures). The predicted [Fe n]/[0 i] ratio assumes 
that the mass in each element of volume was about the same 
between 1998 and 2005, i.e. that the difference in n e between 
the two data sets is caused by an increase in ionization fraction 
between the two epochs. The results are remarkably similar to 
the shock method: an Fe/O ratio lower than solar by a factor 4 
in the IVB interval and a factor 12 in the LVB interval. Again, 
inferred depletions would increase if the total jet density (in- 
stead of x e ) increased between the two epochs. Finally, we note 
that applying an extinction correction to the observed line fluxes 
would decrease the observed [Fe n]/[0 i] ratio, again increasing 
the inferred depletion in both methods. 

To evaluate a possible effect of time variability in the 
[O i] line, we compared t he [O i] fluxes measured in 1998 
by lLavallev-FouQu et et al.l yOOO) against the flux-calibrated 
HST/STIS transverse PV diagram obt ained on 2003 Dece mber 
1 at a distance of -a.' 3 from the star dCoffev et al.ll2007l) . The 
[O i] line fluxes per arcsecond length of the jet in the IVB and 
LVB intervals in the latter data (1.9 x 10~ 13 and 1.3 x 10~ 13 
ergs _1 cm _2 arcsec _I respectively; Coffey, private communica- 
tion) are 80% and 5% higher than their respective values in 
December 1998 at the same distance, averaged over 075 along 
the jet. Using the Coffey et al. [O i] fluxes would thus increase 
depletion in the IVB range at 073 by a factor 2, while depletion 
in the LVB would remain unchanged. However, line fluxes mea- 
sured in the narrow HST slit (of width 071) will also be more 
affected by the precise position of knots at the epoch of observa- 
tion than fluxes averaged over 075 along the jet. Results obtained 
by comparison with Lavalley [O i] data may thus be more reli- 
able despite the longer time span. 

We thus tentatively interpret the low [Fe n]/[0 i] ratios ob- 
served as indication of iron gas-phase depletion by a factor 
=s 3 - 4 in gas faster than -100 km s _1 and by a factor of 10 - 12 
at speeds below -100 km s _1 . This would correspond to about 
33% and 10% of elemental iron in the gas phase, respectively. 
Simultaneous observations of the DG Tau jet involving refrac- 
tory and non-refractory species would be extremely valuable to 
confirm the present findings. The observed trend of higher de- 
pletion at lower flow velocities, also seen in iPodio et al.1 (§009) 
in calcium toward the HH 111 jet, would imply that slower jet 
material in the jet originates in regions farther from the star than 
higher velocity material, or that it has been less shock-processed. 
Implications for launching models are discussed in Section l4~5l 

4.3. Blueshifted jet mass-flux in the [Fe nj emitting 
components 

From our [Fe n] observations we derive estimates of the mass- 
flux rates in the atomic component of the blueshifted jet as a 
function of the distance from the star, and for the two velocity 
intervals HVB and MVB. We used three different methods, the 
first one based on estimates of the jet density and cross-section, 
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the other two based on the [Fe n] luminosity, corrected for deple- 
tion. We detail below the assumptions and parameters for these 
three methods. The resulting mass-fluxes are plotted as a func- 
tion of distance in Figure [9] 

4.3.1 . Method 1 : from jet density and cross-section 

For a jet of cross-section A j and velocity Vj, uniformly filled 
with a density of hydrogen nuclei n#, the mass-flux rate is given 

by 

Mj = fim H n H AjVj, 

where /u, the average weight per hydrogen nucleus, is 1 .4 for a 
cosmic gas with 10% helium. Assuming a circular cross section 
of diameter Dj, the above equation can be re-written as 



100 kms~ 



This method has the advantage of being unaffected by iron 
depletion or by circumstellar extinction. We plot in Figure [9] 
(dashed-dotted line) the mass-loss rates derived from this 
method and computed every 0'.' 1 along the blueshifted jet. 

We used an average deprojected velocity Vj = 300 km s -1 
for the HVB and Vj - 135 km s for the MVB ranges, derived 
from the average radial velocities of V r ^ 220 and 100 k m s _I 
of the spectral components resolved by iPvo et alj ([2003), with 
a jet incl i nation angle with respect to the line of sight of i=42° 
Lavallev (2000). To estimate the HVB cross-section, we took 
D hvb = fwhm Ivb - FWHMl where FWHM HVB is the trans- 
verse [Fe n] emission width measured in the HVB deconvolved 
maps, and FWHM = 0705 - 0708 is the effective resolution 
of these maps. To estimate the cross-section of the MVB flow, 
we subtracted the cross-section of the central HVB component 
nested inside it, i.e. we take D 2 MVB = FWHM 2 MVB - FWHM 2 HVB . 
The effective resolution of the deconvolved maps cancels out in 
this case. 

Estimates of the total density n# along the jet are derived 
from njx e , where the electronic densities n s along the jet are ob- 
tained from our [Fe n] 1.533/1.644 line ratios. For the hydrogen 
io nization fraction jc e , we used linear fits to the values derived 
bv iBacciottil J2002) in 1999 with HST/STIS from optical line 
ratios on similar spatial scales and velocity ranges. One poten- 
tial limitation in using these values may come from time vari- 
abilit y. However, the x e values obtained in 2003 at d = -073 
by Coffey et alJ d2008h agree very well with the 1999 values of 
IBacciottil d2002h . indicating moderate variability on a three year 
timescale in this parameter. Our mass-flux estimates at d — -073 
in the H VB and MVB inter vals differ slightly from those ob- 
tained bv lCoffev et alJ (|2008) using the same method (star sym- 
bols in Figure[9]l, the discrepancy mainly comes on the one hand 
from a jet radius twice as large as our estimate (they do not 
correct for instrumental PSF) and because of the quadratic de- 
pendence, this resul ts in a factor 4. On the other hand, the total 
density obtained bv lCoffev et alJ (2008) is a factor 3 lower than 
ours, which compensates the factor 4 from the jet width. The 
difference in MVB also arises because the authors do not correct 
the MVB width by the area filled by the HVB, which we do. 

4.3.2. Method 2: from [Fe n] luminosity, volume emission 
from uniform jet slab 

The second method derives the mass-loss rate from 
[Fe n] 1.64/im luminosities assuming optically thin vol- 
ume emission from a uniform jet slab. This method, described 



in detail in lHartigan et alJ (1 19951) for optical lines, assumes 
uniform T e , x e and n e and relates the line flux to the emitting 
mass within th e elementary apertur e. Using a 16-level excitation 
model of Fe + (Pes enti et al.ll2003l) . we fitted the emissivity of 
the 1.64/im line as a function of n e for T e 10 4 K as 

/ 3 5 x io 4 \ _I 

j(1.64yum) = 2 x 10~ 17 x 1 + erg s" 1 sr" 1 ion -1 (3) 

\ « e (cm s )j 

(a good approximation for n e ^ 3000 cm -3 ). The mass-loss rate 
is then given by the mass divided by the pixel crossing time 

3.5 x 10 4 \ / Lpan 



M 



1.45 x 10" 8 M o yr _1 (l + 
( (l50 km s =r )( 



n e (cm- 3 )J\lQ- 4 L, 



[Fe/H] 



2xl0 15 cm/ \[Fe/H] G 



(4) 



with /, the size in the plane of sky of the elementary aperture 
along the jet axis, V, the tangenti al flow velocity in th e plane of 
sky, and [Fe/H] = 2.82 x 10" 5 (lAsnlund et alJl2005h . 

We computed the [Fe n] 1 .64/im luminosity in the HVB and 
MVB channel maps in rectangular apertures of projected length 
along the jet axis I, = 0.1" = 2 x 10 14 cm (at 140pc), and width 
±075 across the jet. We adopted the same values of n e as a func- 
tion of distance as in method 1 and the same radial velocities to 
derive V, = V r x tan(/). Finally, we applied the correction for the 
Fe gas-phase depletion estimated in the previous section, using 
the "volume" hypothesis, [Fe/H]/[Fe/H] Q = 1/4 and 1/12 for the 
HVB and MVB respectively. The mass-loss rates derived from 
method 2 are plotted in full line in Figure [9] 



4.3.3. Method 3: from [Fe n] luminosity, emission from shock 
fronts 

The third method to estimate the jet mass-loss rate from the 
[Fe n] luminosity assumes that the line emission in the elemen- 
ta ry aperture is arising f rom shock fronts. It is described in detail 
in lHartigan et al.l (fl995) for the [O i] line. The mass-flux entering 
the shock, M s , is given by 

M s = nm H n H V s A s , 

with id — 1 .4, and V s and A s , the shock velocity and cross-section 
respectively. 

Fro m the planar atomic J-shock models of lHartigan et al.l 
(2004b), we derive a proportional relation between n# V, and the 
emitted [Fe n] 1.64/im line flux per shock area (multiplied by 2 
to account for both sides of the shock) 



' 1.64/an 



9.69 x 10~ 14 x n H V : 



[Fe/H] | _, _ 2 

[F^]7 ,ergs cm 



where we scaled the model fluxes to the recent estimate of so - 
lar abundance [Fe/H] = 2.82 x 10~ 5 dAsplund et al.ll2005h . 
Integrating over the shock area A s yields the following equation 
for the shock mass-flux 



M s = 9 x 10~ 8 M Q yr 



LffW \( [ Fe / H ] 



lO- 4 L /\ [Fe/H] 



(5) 



The jet mass flux is related to the mass-flux entering the 
shock(s) inside the a perture, M s , by the following equation (see 
lHartigan et al][T9^ICabridl2002h 



Mj = M s 



vA( cose \ 



(6) 
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where cos(#) = A j/A s for oblique shocks at an angle 9 with the 
flow normal, and N s h oc k is the number of shock fronts within the 
integration aperture. 

Here we assume perpendicular (cosO = 1) shock fronts lo- 
cated ever y 50 AU along t he jet (m aximum cooling length in the 
models of Hartigan et al.l d2004 b)). This will tend to give an up- 
per limit to the mass-flux. The [Fe n] luminosity in the HVB and 
MVB channel maps is integrated every 0.3" = 50 AU along the 
jet axis, so that N s hock - 1 ■ In addition, following the e stima tes of 
shock velocities derived by lLavallev-Fouquet et al.l d2000h . we 
adopted Vj/V s = 4 in the HVB and Vj/V s = 2 in the MVB 
intervals. Finally, we corrected for the Fe gas-phase depletion 
estimated in the previous section in the shock hypothesis, using 
[Fe/H]/[Fe/H] = 1/3 and 1/10 for the HVB and MVB respec- 
tively. The derived mass-loss rates are represented with a dashed 
line in Figure [9] 

4.3.4. Comparison of the three methods 

We plot in Figure [9] the mass-loss rate estimates in the blue DG 
Tau jet as a function of distance to the star obtained for the three 
methods, with a separate panel for the HVB [-300,-160] km s _I 
and MVB [-160,-50] km s~'velocity intervals. The three meth- 
ods agree reasonably well within a factor 2-5 at all distances 
from the star, lending some confidence to these values, and to 
our depletion estimates. 

For comparison, the mass-fluxes obtained by 
lLavallev-Fouquet et al] d2000l) in their IVB and LVB inter- 
vals using the optical line [O i] and methods 2 and 3 are also 
shown in Figure [9] (gray/green curves). They were scaled up 
by 1.44 to be consistent with the new solar o xygen abundance 
[O/H] = 4.57 x 10~ 4 of lAsplund eTail d2005l) . Because the iron 
depletion was inferred from the [O i] to [Fe n] luminosity ratio 
assuming the same jet density in the two epochs, it is expected 
that both lines should give identical mass-fluxes. The slight 
differences observed stem from the use of slightly different 
velocity intervals (IVB, LVB rather than HVB, MVB), and from 
the fact that the [O i] luminosity is not strictly proportional to 
tihV s in shocks (the slope is steeper than 1), so that method 3 is 
less accurate for [O i] than for [Fe n]. 

Averaging our [Fe n] mass-flux estimates over the central 1" 
gives a mass-loss rate of 1.7 + 0.7 x \0~*M & yr~ x for the HVB 
and 1 .6+0.9X 10"" 8 M o 3?r _1 for the MVB using the three methods 
(the quoted uncertainty range is the 1 -sigma dispersion about the 
mean). Our derived HVB mass-loss rate is c onsistent with the 
lower limit of 1 .2 x \Q^M Q yr l estimated bv lBeck etall d2010l) 
from high-velocity extended HI Bry emission. The total mass- 
loss rate for the [Fe n] emitting atomic gas in the velocity range 
-300 to -50 km s _1 is therefore Mj ^ 3.3 ± 1.1 x lO^M^r -1 , 
with approximately equal contribution from the HVB and the 
MVB. This value is a lower limit to the total mass-loss rate of the 
atomic component because the [Fe n] emission does not probe 
the whole range of velocities seen at optical wavelengths. 

4.4. Accretion rate and ejection-to-accretion ratio in DG Tau 

iMuzerolle et al.l (fl998) have shown that the luminosity of the 
Bry line in a large sample of T Tauri stars is well correlated 
with the accretion luminosity onto the star inferred from mod- 
eling of the optical-UV excess. During our SINFONI run, we 
also obtained K-band data including the Bry line. We find that 
any spatially Bry extended component contributes at most 2 % 
of the total line flux integrated over all velocities. This is consis- 
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Fig. 9. Values of the mass-loss rate for the atomic blueshifted jet 
as a function of projected distance to the star and for the two ve- 
locity components: HVB (top panel) and MVB (bottom panel). 
The values derived from the [Fe n] 1.64/vm observations using 
the three methods discussed in the text are plotted with black 
curves: jet density and cross-section (dashed-dotted curve), uni- 
form slab ( full curve), shock front (dashe d curve). The values 
obtained by Lavall ev^Fouquet et al.l (12000) using the optical line 
[O i] and methods 2 and 3 are shown for compariso n (gray/green 
curve s). Star symbols show the values obtained bv lCoffev et al.l 
(120081) from an optical HST/STIS transverse slit at -a/3. 



tent with the recent finding o f lBeck etal.l d2010). who show that 
although extended high-velocity Bry emission is detected along 
the DG Tau blueshifted jet, the extended component makes up 2 
% of the total (velocity integrated) line flux. We therefore make 
the assumption below that all the HI Bry flux in DG Tau is re- 
lated to the accretion process. From our flux-calibrated data we 
derived a K-band continuum of 7.1 mag and a Bry equivalent 
width of 5.6 A corresponding to a line brightness of 3.9 x 10~ 16 
W rrT 2 and line luminosit y of 2.5 X 10~ 4 L n afte r extinction cor- 
rection by Ay = 1.6mag (|Gullbring et al.l r2000). Using the em- 
pirical fit of IMuzerolle et al.l(ll998l) . we derived an accretion lu- 
minosity L acc = O.77L0 for DG Tau in October 2005. Assuming 
that the luminosity arises from an accretion shock at the base of 
an accretion column anchored at the inner disk truncation radius 
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Ri, the mass accretion rate onto the star is then obtained from 
dGullbringetal.ll 19981) 



~gm7 



(7) 



Taking = 5/?* as in lGullbring etal.1 (H998), and adopting 
for DG Tau a stellar radius of R+ =23 Rq and a stellar mass of 
M*=0.7 M^ dHartigan et alJll995h we derived an accretion rate 
of M acc = IP' 7 Mpyr" 1 . With a similar Bry flux (4.5 x 1(T 16 W 
Beck et all dfo 1 0) derive M acc 



m 



9.6 x 1CT S Moyr -1 using 
Ay=3.3 mag and a stellar mass of 2.2 M for DG Tau from 
IWhite & Hillenbrand! d2004l) . We note that this value of M acc 
only probes the instantaneous accretion rate from the inner disc 
onto the stellar surface in October 2005. It is not necessarily rep- 
resentative of the mean accretion level pertaining to the period 
2002-2004 when the jet material probed by our SINFONI obser- 
vations was ejected (see Fig.[7]l. 

Indeed, diagnostics of the accretion shock, such as the Bry 
and Ca II line strength or the optical veiling, are known to be 
highly variable in T Tauri stars on sca les from days to years , 
especially in stars as active as DG Tau. iMuzerolle et all (1998) 
measured in 1998 a Bry equivalent width 3 times larger than 
ours, EW(Bry)=14.1 A, and obtained a 3 times higher accretion 
rate, assuming a typical K-band magnitude for DG Tau K= 6.74 
(avera g e of 17 measurements, provided in dKenvon & Hartmannl 
ll995l) V lGullbring etal.1 d2000h derived M acc = 5 x 10~ 7 M yr -1 , 
by fitting an accretion shock model to the Balmer and UV excess 
continuum emission measured in 1996. Finally, the mean corre- 
lation between Call line luminosity and accretion rate in T Tauri 
stars suggests M acc ^ 2 x 10" 7 M yr _1 in DG Tau in 2002-2003 
dMohantv et al]l2005l) . 



Below we therefore adopt M a , 



(3 ± 2) X 10-' Moyr 1 as 



the possible range of accretion rate prevailing in DG Tau. This 
is sig nificantly lower tha n the M acc = 2 x 10~ 6 M yr~' obtained 
bv lHartiganetaild 19951) from an analysis of the optical veiling 
in 1988-1989. However, the method of HEG95 systematically 
tends to overestimate acc retion rates in T Tauri stars by an order 
of magnitude on average dGullbringetal.1119981) . 

From the range of mass-flux rates derived in the [Fe n] jet 
(see Section l4~3l l. we thus inferred for the blueshifted jet lobe 
a ratio of ejection to accretion rates of 0.04-0.4 for the [Fe n] 
emitting flow, with about equal contributions from the MVB 
and HVB ranges. Most of the uncertainty in this ratio (a fac- 
tor 5) arises from the mass-accretion ra te. The ratio is substan- 
tially h igher than previously reported bv lLavaliev -FouQuet et alJ 
(2000), m ainly because they ado pted a 10 times higher accretion 
rate from lHartigan et ail d 19951) . It is compatible with the aver- 
age ejection to accretion ratio for one jet lobe of ^ 0.1 in atomic 
T Tauri jets obtained by ICabritl d2007bl) when c omparing [O i] 



mass- flux rates from spatially unresolved data (Hartigan et al 
1995) with revised mass-accretion rates from iGullbring et al 
( 1998). Our mass-flux estimates therefore do not appear to suffer 
from a strong systematic bias. 



4.5. Implications for jet launching models 

Our results on mass-fluxes, kinematics, and iron depletion in the 
DG Tau jet have several implications for proposed jet origins. 

First, we note that the observed velocities in excess of -50 
km s in the blueshifted [Fe n] flow are much too high to be pro- 
duced by disk photo-evaporation. Indeed, in photo-evaporated 
flows terminal velocities reach at most 3-4 x the sound speed 



cs at the base of the flow dFont et a l. 2004). Maximum tempera- 
tures of =i 10 4 K reac hed at the inner disk surface through EUV 
and X -ray irradiation (lErcolano et al.l l2009: Gort i & Hollen bach 
2009) would thus lead to terminal velocities < 30-40 km s , 
much lower than observed in [Fe n]. 

One might wonder if the MVB flow could trace matter en- 
trained and accelerated in a mixing layer between the high- 
velocity jet (traced by the HVB) and the outer photo-evaporated 
disk wind. The most recent computations of disk p hoto- 
evaporation dOwen et al.ll2010t iGorti & Hollenbachll2009l) . in- 
cluding X-ray irradiation, indeed predict mass-flux rates that 
could reach 10 9 -10 8 M yr~' over the entire disk, sufficient to 
supply the MVB flux. However, the same models show that most 
of the photoevaporated mass-flux arises from disk radii r > 5 AU 



(IGorti & Hollen bach 



2009UOwenetal ll2010l) . The correspond- 
ing streamlines (see Ow en et al. Fig. 3) have cylindrical 
radii > 20 AU for altitudes above the disc z > 30 AU, signif- 
icantly wider than the derived HVB jet radius (see Figure [3); 
therefore, one expects very little or no entrainment of the photo- 
evaporated flow by the HVB jet. 

Second, we note that the gas-phase depletion of iron, if con- 
firmed, rules out an origin in the stellar surface, at least for 
the [Fe n] MVB emitting material. Indeed, recent determina- 
tion of photospheric [Fe/H] f or young stars in Taurus are com- 
patibl e with solar abundance dSantos et al.ll2008l iD'Orazi et al.l 
1201 ll) . Lower [Fe/H] ratios (^ 0.2 solar) have been derived 
in X-rays tracing coronal material but in these observations 
the [Fe/O] abundan ce ratio is 0.4-0.8 solar dScelsiet al .1120071: 
Giint her et al.l l2007h . These abundance patterns cannot repro- 
duce the low [Fe/O] ratios (^ 0.1 solar) inferred in the MVB 
emitting gas, indicating that the observed Fe gas phase deple- 
tion is likely caused by depletion onto dust grains. Our observed 
Fe depletion pattern also excludes the fact that the medium- 
velocity conical sheath surrounding the jet could trace a cocoon 
of matter ejected sideways from the narrow high-velocity beam. 
Sideways ejected material should be equally or less depleted 
onto dust (more shock-processed) than the faster central beam, 
whereas we observe the opposite, with a 3 times larger deple- 
tion at lower velocity (see Section l4~2l ). It is difficult to pinpoint 
more precisely the launch radius from the depletion alone; the 
dust sublimation radius in DG Tau appears to b e r su b - 0.14 
AU f rom K-band interferometric measurements dAkeson et al.l 
120051) . However, the K-band size may be o nly an upper limi t 
to r su b owing to the effect of scattered light dPinte et al.l l2008). 
Ejection from near corotation (r corot = 0.06 AU in DG Tau with 
a stellar period of 6.2 days and assuming a O.7M star) is thus 
not totally excluded at this stage, pending more detailed mod- 
els. Ejection from slightly beyond 0.1 AU is not ruled out ei- 
ther, because sublimation starts further out for smaller grains, 
and only 10%-33% of the dust appears to have been destroyed 
in the [Fe n] jet. In addition, r su b could be larger than the K-band 
radius if the inner disk of DG Tau is optically thick in the gas 
continuum. 

Below, we confront the DG Tau jet properties with pre- 
dictions for proposed ejection processes, including the X-wind 
model and an extended MHD disk- wind model. 



4.5.1. The X-wind model 

The X-wind model assumes that the disk is not intrinsically mag- 
netized and is threaded only by the stellar dipole. It posits that 
stellar fieldlines are bunched up at the corotation radius and are 
partly opened in a fan-like configuration along which a centrifu- 
gal wind, carrying off the angular momentum from the funnel 
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flow, is assumed to be launched. The semi-analytical wind so- 
lution is characterized by an axial density enhancement (the jet) 
surrounded by a wide-angle flow, where the spe ed is nearly con- 
stant across all streamlines (Sha ne et al.l[l9 98). This predicted 
constant flow speed stems from the cylindrical shape of the 
Alfven s urface, which a ppears to be a generic property of the 
X-wind dCai et al] 12008 ). Because all streamlines are launched 
from the same radius ro = r oorot , their magnetic lever arm param- 
eter A = (r^/ro) 2 is also the same, and so is the terminal speed af - 
ter magneto-centrifugal acceleration dBlandford & Pa vnal 19821) 



Vp = V2/1-3V, 



Kep 



(ro). 



(8) 



The steep velocity decrease across the DG Tau blueshifted jet 
seen in Figs 4 & 5, from V=-200 km s~'in the HVB for angles 
to the jet axis less than 4° to V=-100 km s^'in the MVB for an- 
gles to the jet axis of 14° and down to -30,-40 km s" 1 at w ider 
angles in [O i] and [S n] dBacciottfeOOl ICoffev et al.ll2007l) . is 
clearly not compatible with this characteristic property. The X- 
wind model could account for the [Fe n] HVB component only, 
but this would require that the MVB flow originates in an MHD 
disk wind (because it cannot be explain ed by photo-evapor ation, 
see above). However, as discussed in iFerreira et alJ (120061) . it is 
very unlikely that an X-wind and an extended MHD disc-wind 
will coexist because this would imply a hole in the disc mag- 
netic flux distribution between these two components. Therefore 
the X-wind model as currently envisioned does not explain the 
kinematics in the [Fe n] DG Tau jet. 

4.5.2. Conical wind from the disk-magnetosphere boundary 

Recent numerical simulations of the interaction between a stel- 
lar dipole and a disk with no net magnetic field yield a situation 
somewhat different fro m the (semi-analytical) X-wind model 
( Rom ano va et al . 2009). As the stellar fieldlines are bunched up 
by accretion and stretched, a narrow conical wind is formed 
along the neutral line (rather than in the angular sector below 
it), and powered by the magnetic pressure force (rather than by 
the centrifugal force as in the X-wind). In the standard (non- 
propeller) regime relevant to DG Tau, where the star rotates 
slowly, the total mass-flux rate in the conical wind is about 10%- 
30% of the accretion rate onto the star, compatible with the ob- 
servational constraint. A strong transverse velocity gradient is 
also present, with faster matter from the stellar magnetosphere 
on the inside and slow matter from the disk on the outside. 
However, most of the wind mass reaches a low poloidal speed 
of about 40 km s _1 for the CTTS dimensional scalings appropri- 
ate in DG Tau: v„ ^ .2vq with Vo = 197km s in Table 1 of 
iRomanova elaTl d2009l) . And the half-opening angle of 30°-40° 
is much wider than the 4-15° observed in [Fe n], Therefore the 
conical wind as modeled in this work appears both too slow and 
too broad to explain the [Fe n] jet observations, although it could 
match the LV sh eath surrounding the jet observed in optical lines 
(lBacciottill2002l) . In this framework, the whole [Fe n] jet (HVB 
and MVB) would have to originate mainly in the stellar surface. 
This would contradict with the observed iron gas-phase deple- 
tion in the DG Tau jet, indicative of a substantial dust content. 

4.5.3. Extended MHD disk wind 

Protostellar disks may indeed harbor a substantial magnetic flux, 
advected during the initial infall phase of star formation. If 
the field is slightly below equipartition, rotation remains quasi- 
Keplerian in the disk midplane and angular momentum can be 



efficiently extracted by magnetic torques to launch a powerful 
wind from the disk surface, as shown by analyti cal and numer- 
ical 2D solutions of the vertical disk structure dFerreiral ll997; 
IZanni et al.ll2007l) . iPesenti et al] d2004h found that a self-similar , 
steady MHD wind from a warm disk dCasse & Ferreiral 2000). 
with a lever arm parameter A = (rA/ro) 2 - 13 could fit the tenta- 
tive jet rotation signatures observed in the low-velocity [O i] flow 
(Bacciotti 2002). The same type of m odel can also reproduce 
observed jet widths in T Tauri stars dRav et al.l 12007). We now 
check whether this model could also explain the speed, mass- 
flux, and depletion in each of the [Fe n] velocity components. 

We may use the range of flow speeds in each [Fe n] ve- 
locity component to constrain the range of ejection radii in the 
considered disk wind model. The asymptotic velocity given in 
Eq. [8] is reached only fa r from the star in self-similar solu- 
tions dFerreira et alj|200q) . We therefore use here the exact pre- 
dicted flow speeds at an altitude of 50 AU above the disk for 
the A — 13 model, plotted as a thick black curve in Fig. 3 of 
IFerreira et alJ (2006) (note that velocities in that graph are scaled 
by 1/ yfWjM^). The radi al velocity range of 160 to 250 km s -1 
for the HVB component dPvo et al.l 120031) corresponds to de- 
projected flow speeds of 210-350 km s'and to launch radii of 
0.08 AU to 0.25 AU for a stellar mass M* = 0.7M Q . The up- 
per value is consistent with the outer launch radius of 0.2-0.5 
AU for the HVB gas deri ved from tentative rotation signatures 
in optical/near-UV lines dCoffev et al.ll2007h . Similarly, the ra- 
dial velocity range of [-50,-160] km s _1 of the MVB, and de- 
projected speeds of 70-210 km s _1 , would imply launch radii of 
0.25-1.5 AU. The observed Fe depletions could be explained if 
10% of dust is already destroyed around 1 AU and 33% around 
0.2 AU. Detailed modeling including a dust size distribution and 
opacity of the gas disk are needed to se e whether this is c ompat- 
ible with K-band size measurements ( Akeson et al. 20051). 

The ejection-to-accretion ratio for one jet lobe in a self- 
similar MHD disk wind powe red by angular mom entum extrac- 
tion from the disk is given by dFerreira et al. 2003) 



Mj 



M acc 4(A- 1) 



log(r e /r ; ) = 



1 



2(A-l) 



log(V ma x/V min ), 



(9) 



where r e and r, are the external and internal radii of the wind 
launch zone, and V m i n and V max the corresponding minimum and 
maximum flow speeds. With A — 13, and values of V m ; n and V max 
as above, one expects an ejection/accretion ratio for one jet lobe 
of 0.04 in the MVB and 0.02 in the HVB intervals, respectively. 
This is compatible within uncertainties with the mass-fluxes esti- 
mated here (see Section l4~3l if the DG Tau accretion rate during 
the period (2002-2004) was ^ 5 x 10~ 7 M yr _1 , near the up- 
per end of the range reported in recent works (see Section l4~4l >. 
Clearly, a closer monitoring of accretion and ejection signatures 
in DG Tau over a few years is required to estimate more accu- 
rately the ejection/accretion ratio in the DG Tau jet, and provide 
a more definitive test. 

We stress that the above comparison must be taken as in- 
dicative only. First, we have discussed the implications in the 
context of a specific disk wind solution with A-\3>. A smaller 
A solution may be more compatible with the mass flux con- 
straints. A thorough investigation of possible disk wind solu- 
tions compatible with all observational constraints should be 
conducted. Secondly, the above comparison assumes that the 
periodic knots/bubbles do not significantly modify the overall 
steady-state disk wind kinematics and mass-flux (even though 
they will locally modify its ionization level through shocks). The 
period of knot generation of 2.5 years in DG Tau is about 140 
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times the star rotation period of 6.2 days, i.e. much longer than 
the orbital timescale in the inner disk, but comparable to the or- 
bital period at 1 .5 AU. Therefore, the steady-state disk wind as- 
sumption could be valid for the HVB, but only marginal for the 
MVB flow. The formation of large bubbles in the wake of knots 
may further perturb the structure of the outer disk wind. 

Large-scale flares over similar long timescales, followed by 
a re-arrangement of outflow density and velocity distribution, 
were recently observed in long-term numerical simulations of an 
extended dis k wind interacting with the stellar magnetosphere 
dFendtll2009l) . This scenario could be promising to explain the 
DG Tau jet properties, although simulations over larger spa- 
tial scales would be needed to compare with our observations. 
Alternatively, a stellar magnetic cycle could be at the origin of 
the time variability observed in the launching process or the pe- 
riod of 2.5 years might correspond to an excentric binary com- 
panion perturbing the inner disk at periastron passage. Studies of 
jet wiggling and interferometric imaging of DG Tau could help 
in testing this last alternative. 

4.5.4. Episodic magnetic bubbles 

Based on re cent laboratory exp eriments of radiatively cooled 
plasma jets, ICiardi et all d2009t) propose yet another possible 
ejection mechanism for the DG Tau jet. Starting from a highly 
wound-up field configuration with \B$\ oc r _I » B p and a strong 
axial current, two outflow components are generated: a magnetic 
bubble (or cavity) accelerated by gradients of the magnetic pres- 
sure and surrounded by a shell of swept-up ambient material, and 
a magnetically confined narrow jet on the bubble axis. This pro- 
cess is akin to th e magnetic tower modeled numerically e.g. by 
iKatoetalJ (12004 . although it is followed over longer timescales 
and with higher resolution. Repetitive outflow cavities are pro- 
duced as the plasma refills in the launch region. The new cavities 
are confined by the matter left over from previous episodes, so 
that collimation properties gradually become independent of ini- 
tial conditions. The morphology is suggestive of that of DG Tau, 
with a collimation angle of < 10° for the high-density clumps 
along the axis. However, the bubbles already inflate very close 
to the disk so that the conical geometry of the MVB flow is not 
well reproduced. Additional confinement of the bubbles by a sur- 
rounding disk wind may be required to explain this property. 
Further work in this direction would be extremely valuable, as 
would numerical estimates of the required \B$\ in the DG Tau 
disk. 



5. Conclusions 

- We observed an onion-like kinematic structure in the 
blueshifted jet with a conical medium-velocity (MVB) flow 
at V a -100 km s _1 surrounding a more tightly collimated 
high-velocity (HVB) beam at V a -200 km s _1 . The velocity 
fluctuations and electronic density distribution along the jet 
show little variation between 1998 and 2005. However, little 
[Fe n] emission is observed from the low-velocity (LV) flow 
component at a -30, -40 km s _1 reported in optical lines at 
30 AU from the jet axis. 

- By comparing the intensities in [Fe n] and [O i] in the 
DG Tau jet at two epochs, we infer evidence for iron de- 
pletion by a factor 10 at speeds below 100 km s , and 
by a factor 3 in the high-velocity range. Simultaneous data 
are required to confirm this result. Mass-flux rates for the 
blueshifted jet lobe inferred from three independent meth- 



ods are 1.6(±0.8x lO^Moyr" 1 in each of the MVB and HVB 
ranges, representing 0.02-0.2 of the disk accretion rate each. 

- The counterjet is dominated by a striking bubble struc- 
ture. A similar faint bubble structure is also detected in the 
blueshifted lobe. The radial and tangential velocity in the 
red lobe is lower than in the blue lobe by a factor a 0.7. 
The observed central obscuration, if caused by extinction by 
the disk, favors the similarity solution f or the DG T a u disk 
surface density distributions derived by llsella et ail d2010h 
rather than a sharply truncated power law. 

- The two new knots present in the blueshifted lobe indicate 
the presence of a variability timescale of 2.5 year s over the 
last 10 years, two times shorter than estimated by iPvo et all 
d2003l) . and 140 times longer than the stellar rotation pe- 
riod. The current bubble structures appear to ha ve formed 
in the wake of the high-velocity knot observed bv lPvo et al.l 
(2003), ejected in 1999. Observations more closely spaced 
in time would help to clarify the knot/bubble dynamics and 
the origin of these quasi periodic eruptions. Searches for pe- 
riodic changes in mass accretion rate and for jet wiggling 
would also constrain the process responsible for these tran- 
sient events. 

- The velocities of the [Fe n] MVB flow are too high to 
be accounted for by disk photo-evaporation models. Origin 
of the MVB component in entrainment of an outer photo- 
evaporated disk wind seems also unlikely owing to the de- 
rived narrow HVB jet widths. The strong iron depletion 
also rules out an origin in a stellar wind, or in sideways 
ejected matter from the (3 times less depleted) HVB beam. 
In addition, neither the classical X-wind model (Shang et al. 
1998; ICai et al.l 120081) nor the conical wind found in re- 
cent simulations of ste llar magnetosphere-disk interaction 
dRomanova et alj[2009l) are able to reproduce the kinemat- 
ics and opening angles in the [Fe n] MVB. A dusty MHD 
disk wind appears to be the most plausible origin for this 
component. 

- The extended MHD disk-wind solu tion proposed to fi t opti- 
cal rotation signatures in DG Tau ( iPesenti et alj|2004l) com- 
pares more favorably with observations. It predicts launch 
radii of 0.08-0.25 AU and 0.25-1.5 AU, respectively, for the 
HVB and MVB velocity ranges, consistent with the higher 
iron depletion at lower speed. The predicted mass-fluxes also 
agree with observed values if the accretion rate was rela- 
tively high in DG Tau over 2002-2004, ^5x 10' 1 M e yr~ l . 
However, contribution from magnetospheric reconnexions or 
magnetic towers may not be negligible, modifying the pre- 
dicted flow structure. Spectro-imaging monitoring of accre- 
tion and ejection signatures in DG Tau over a few years are 
essential to elucidate the ejection processes at work in this 
prototypical source. 
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